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In order to determine the origin of discontinuities which arise when the semiclas¬ 
sical propagator is employed to describe an infinitely long and infinitesimally thin 
solenoid carrying magnetic flux, we give a systematic derivation of the semiclassical 
limit of the motion of an otherwise free charged particle. Our limit establishes the 
connection of the quantum mechanical canonical angular momentum to its classical 
counterpart. Moreover, we show how a picture of Aharonov-Bohm interference of 
two half-waves acquiring Dirac’s magnetic phase when passing on either side of the 
solenoid emerges from the quantum propagator, and that the typical scale of the 
resulting interference pattern is fully determined by the ratio of the angular part of 
Hamilton’s principal function to Planck’s constant. The semiclassical propagator is 
recovered in the limit when this ratio diverges. We discuss the relation of our results 
to the whirling-wave representation of the exact propagator. 
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I. INTRODUCTION 


In a series of thought experiments, Aharonov and Bohrn demonstrated that charged 
particles are influenced by electromagnetic potentials, even if their wave function vanishes 
wherever the associated electromagnetic held is nonzero [Tj. A controversial debate on 


whether such effects existed has since been resolved by their experimental verification | 2 J. 


Phase shifts due to the presence of inaccessible magnetic flux have been confirmed in electron 
microscope experiments as well as by means of electron holography [ 6 ], and with the 
same technique where the flux was additionally shielded by a superconducting layer, to 
further minimize the influence of stray magnetic fields [?]. The simulation of Aharonov- 


Bohrn phases in cold atom systems trapped in optical lattices [ 8 ] is now experimentally 


feasible [9], and the effect plays an important role in mesoscopic electron-optical systems pTO]- 
[ 12 ], where an interesting manifestation for charged identical particles [Till fl 6 ] has recently 
been confirmed in an electronic Hanbury Brown and Twiss interferometer DU- 

An elementary realization of the effect, first put forward by Ehrenberg and Siday [IB], 
which also appears in the publication by Aharonov and Bohm |T], is now described in 
several textbooks (see, e.g., [IHJEO]). The thought experiment relies on a semiclassically 
inspired description in which the particle wave reaches a detector via mutually exclusive 
classical alternatives, 7 / and 7 n in figure [lj whereupon each contribution, due to a classically 
inaccessible magnetic flux $, acquires Dirac’s magnetic phase factor |2Tj. The probability 
amplitude which results for the particle wave emanating from the source i to reach a specific 
point on the screen s is then given by 



( 1 ) 


where K^ I/U denotes the contribution from the first/second path in the field-free case, and 
A is the vector potential giving rise to the magnetic flux d>. Taking Dirac’s magnetic phase 
acquired along the first path out of the bracket in the second line in (JT|) corresponds to a 
reversal of this path in figure [I] Then the vector potential which comes with the second 
term in the second line in ([Tj) is integrated along a contour enclosing the solenoid such 


that invoking Stoke’s theorem produces a phase proportional to the magnetic flux Thus 
the interference fringes on the screen are periodically shifted as a function of $, while the 
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FIG. 1 . Elementary manifestation EEBH2D1 of the Aharonov-Bohm effect in a double-slit thought 
experiment in which charged particles, emanating from the source i to reach the screen s, pass either 
side of a thin solenoid placed behind the slit-system, between the two openings. The probability 
amplitudes associated with the two topologically distinct paths 7 j and 7// acquire Dirac’s magnetic 
phase factor, giving rise to a phase difference which amounts to the flux enclosed by the curve 
which ensues after a reversal of either path. 


envelope of the pattern remains invariant, which is confirmed by experiment M. None of 
the sources [Tl fTHl - EU] of the above thought experiment, however, state how such a description 
is related to Schrodinger theory. 

A precise relation assigning Dirac’s magnetic phase to classical trajectories is estab¬ 
lished by the semiclassical approximation of the quantum propagator by Van Vleck and 
Gutzwiller [221122J 


/C (T", x! \ t", t') = 


y/2mh g cl 


det 


d 2 TZ(x", x'\ t", t') 
dx'dx" 


x exp ( jlZ(x", x ; t", t') — i^L— ), 


( 2 ) 


which constitutes an approximation of the probability amplitude of a particle starting out 
at x(t') = x' to be later found at x(t") = x". Here / denotes the number of degrees 
of freedom, and the sum runs over all classical trajectories x c i which connect x' and x", 












4 




(a) Process enclosing the flux 


(b) Displacement of the endpoint 


FIG. 2 . Charged particle propagating from x' to x" in the presence of a magnetic string carrying 
magnetic flux <f>, and in the vicinity of an infinitely extended wall, (a) As long the two possible 
trajectories, 7 / unperturbed, and 7 u once reflected at the wall, enclose the magnetic string, the 
semiclassical propagator 0 generates a probability amplitude similar to @ as encountered in the 
elementary manifestation of the Aharonov-Bohm effect, (b) If the endpoint x" is displaced such 
that the closed curve ensuing after the reversal of either path no longer contains the magnetic 
string, the amplitude (jT]) abruptly changes into ([ 8 ]) which is no longer a function of the flux $. 
This renders (jT]) discontinuous whenever a classical trajectory passes directly through the magnetic 
string. 

entering Hamilton’s principal function 


lZ(x", x'\ t",1') = / £ ( x(t), x(t), t ) dt, 


(3) 


where £ is the Lagrangian of the system (for the origin and significance of the additional 
phase n cf. reference [23]). Thus each classical trajectory in (j2j) contributes a phase de¬ 
termined by Hamilton’s principle function 1Z in units of Planck’s constant. If £ contains 
a vector potential A, the contribution of the latter constitutes precisely Dirac’s magnetic 
phase factor. 

In their original publication Aharonov and Bohrn introduced the vector potential |1] 

r $ 

-4 A r C r -|“ Am&in “ 

Zrcr 


( 4 ) 
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of an infinitely long and infinitesimally thin solenoid, now called a magnetic string, which 
gives rise to a magnetic field confined along the 2-axis. We restrict the discussion to the 
x-y-plane in which the magnetic string thus pierces the origin. The vector potential Q 
generates the second term in the Lagrangian 


f-'x X^Pi 


(5) 


wherein Cq contains all remaining terms when no current runs in the solenoid, and where 
we have introduced the abbreviation 


X = ~ 


e $ 
c 2n 


( 6 ) 


The second term in (J5]) does not contribute to the Euler-Lagrange equations of motion, 
reflecting the fact that classically the particle does not feel any Lorentz force. Yet this term 
contributes to Hamilton’s Principal function 


n x = n 0 - x 


dip 


(7) 


a magnetic phase factor which is proportional to the angle swept by the respective classical 
trajectory [23] - 

If the physical process is such that there exist exactly two trajectories starting at x' and 
ending at x", which enclose the magnetic string, insertion of ([7]) into (j2j) directly generates 
a probability amplitude similar to ([Tj) . A conceivable situation which gives rise to such a 


process is depicted in figure [2a] where in addition to the flux the particle propagates in the 
presence of a wall that extends infinitely in the two-dimensional plane. Another possibil¬ 
ity [25]) where the classical trajectories which enclose the flux are furthermore differentiable, 
arises if the particle in addition to the magnetic string propagates under the influence of 
a Coulomb attractor Z. This also gives rise to at least two classical trajectories 1251 , as 
depicted in figure [3aJ Common to both situations is that the associated semiclassical propa¬ 
gator (J2| is a discontinuous function of the dynamical variables x' and x". For the propagator 
associated with the wall, this can be seen by slightly displacing the endpoint of the propa¬ 
gation process such that the closed curve which ensues by following —7/ after 7 u no longer 


contains the flux, as depicted in figure [2b] For the particle in the Coulomb potential the 
same follows after a rotation of initial and endpoint of the process about the center of at¬ 
traction as depicted in figure [3bJ Thus the probability amplitude abruptly changes from ([!]) 
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(a) Process enclosing the flux 



FIG. 3. Charged particle propagating from x' to x", in the presence of an (infinitely massive) 
Coulomb attractor Z, as well as of a magnetic string carrying magnetic flux <&. The Coulomb 
potential gives rise to at least two trajectories [25] from x! to x ", 7 / connecting both points in a 
rather direct way, and 7 // passing in a sharp turn behind the center of attraction. The classical 
trajectories are invariant under rotations of x! and x" about Z such that the process in (a, adapted 
from [ 22 ]), enclosing the flux T, can be transformed into the process depicted in (b), in which the 
trajectories no longer run around the magnetic string. Thus also in this scenario the semiclassical 
propagator © features discontinuities whenever a classical trajectory passes directly through the 
flux. 


to 


K = exp —- 

he 


1 e 


'11 


A-dx (A'J' + K~<") 


( 8 ) 


whenever a classical trajectory runs directly through the magnetic string. 

The principal intention of the present article is to elucidate the nature of such disconti¬ 
nuities which arise when the semiclassical propagator is employed to describe inaccessible 
magnetic flux. To this end, we will introduce a novel semiclassical limit to determine the 
semiclassical asymptotics of a particle that propagates otherwise freely in the presence of a 
magnetic string. 

The paper is structured as follows: In section [IT] we present established exact representa¬ 
tions of the Aharonov-Bohm propagator, from which we derive our novel semiclassical limit 


in section III Section IV contains a discussion of the interference process implied by the 
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asymptotic limit. In section [V] we elaborate on the relation of onr asymptotic expressions 
to the so called whirling-wave representation of the exact propagator, before we conclude in 
section ED 


II. EXACT SOLUTIONS AND REPRESENTATIONS OF THE 
AHARONOV-BOHM PROPAGATOR 


(9) 


If no further potentials act on our charged particle except for a magnetic string, the 
Hamilton operator associated with the vector potential Q reads 

1 N2 

x 2m y<9r 2 r dr r 2 
which has the eigenfunctions 

i>u (r, (p) = (-i) ]l+x/hl -4 =J\i+ x/ n\ ( kr ) e%Up > 1 e z , k e M+, (10) 

V 2tt 

where J denote the Bessel functions of the first kind m From ( ]10| ) Aharonov and Bohm 
construct their wave function [T] 


OO 

V’x {r, <p) = J\i+ x /h\ (kr) e ll(p , (11) 

l= — OO 

which describes the scattering of a plane wave incident on the magnetic string from the 
positive ^-direction. The eigenfunctions (10) give rise to the completeness relation [27] 

r / ff f\ -j OO /*00 

hirr b W" -v r ) = i r '£, dkk (fcr " ) J l'+*/«! ( kr ‘) e “( 12 ) 

Vr"r' z — y 0 

With the help of the latter one obtains the quantum propagator 

: it J o dkkex P {t " ~ *')) J iuxMi ( fcr ") J i 

( //2 i /2 \ 00 

rmr + r \ 

h~2 t"-t’ J 2s \ l+ xM 

' l=—o o 


AT = 


i mr"r’ 
"ht" -t' 


V+x/*l (U') 


(13) 


where / denote the modified Bessel functions of the first kind m The Bessel functions of 
the first kind J and the modified Bessel functions of the first kind / are related by 


(~i) m J m (x) = I m (-ix), 


( 14 ) 











such that, after the identifications 

p ~ ip" — p\ hr ~ z := 


m r"r' 


h t" — t' ’ 


(15) 


and apart from the overall factor in the second line in (13), the mathematical structure of the 


Aharonov-Bohm wave function ( [II] ) and of the Aharonov-Bohm propagator (13) coincide. 
Therefore all following considerations apply to both objects interchangeably. There exists 


an integral representation of (13), which reads 


K x = K () [ exp ( -~x (<p" ~ p') 


1 

-( 

2 




sin ((^X+XA) 7T) 

and which will serve for later numerical comparisons. Here iA 1 ) denote the Hankel functions 


poo 

x / e" 


of the hrst kind ra. and for the validity of ( |16| ) the angular difference in the hrst term must 
be constrained to values \p" — <p'\ < i r. The integer l x is determined by 


o <i x + x/h < 1, 


and Ao denotes the quantum propagator of the free particle 


(17) 


AV, = 


m 


i m (x" — x') 2 


2itih (t" — t') ^ l h 2 t" — t' 


(18) 


For half-integer flux the propagator ( I3fl6 ) can be written in closed form [lj, which 
we state here as 


dx exp ( ix 2 ) 


K x = K 0 ( exp (~pc W ~ </?')) ~j= [ 

\ \ ^ J w ITT J —y/2z COS ((ip" — if')/ 2) 

f i \ l r-V2zcos((<p"-<p')/2) 

+ exp ( - JX W" - A T 27r) ) -J= J dx exp (ix 2 ) ), 


(19) 


where the sign in the exponential prefactor of the second line is given by that of p" — p'. 


III. SEMICLASSICAL APPROXIMATION AND THE ASYMPTOTIC LIMIT 


Since the magnetic string does not influence the classical motion of the otherwise free 
particle, only one straight line trajectory classically connects x' and x" in the particle-flux 
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line system described by ([9]). Hamilton’s principal function ((7|) of this classical trajectory, 
which we for later purposes express in polar coordinates, reads 


m r" 2 + r' 2 — 2 r"r' cos (p" — p') 
~2 ““ 


X 


dp, 


( 20 ) 


where the first term is the free-particle contribution, and the second term represents the 
contribution of the magnetic string giving rise to the vector potential Q. The canonical 
angular momentum, which is a conserved quantity along the classical trajectory, will play 
an important role in our following analysis and is given by 

,sin (ip" — p') 


pip = mr r 


t" -1' 


X■ 


( 21 ) 


Insertion of (20) into (J2| produces the semiclassical approximation 


K x = /Co exp 


h 


X 


dp 


( 22 ) 


where the semiclassical propagator of the free particle /Co coincides with the exact expres¬ 


sion (18). Due to the restriction imposed on the angular difference which appears in the first 


term of (16), this first term coincides with (22). The latter features a discontinuity whenever 
the classical trajectory associated with the propagation process directly passes through the 


magnetic string (in (16) this discontinuity is compensated by the second term TO, exactly 
as in our example systems of section [TJ Our main purpose here is to elucidate the origin of 
these discontinuities. To this end, we introduce a novel semiclassical limit which transforms 


the full solution (13) for the otherwise free particle into the semiclassical approximation ( 22 ). 
The procedure [26] is based on the exact asymptotic expansion 


Iu(x) 


(- l) k Y{v + k + \) 
Tt c^ {2x) k k\T(v-k + \) 


E 


exp (—x — (y + |) 7 n) EE 1 T(u + k + \) 


V2- 


TTX 


E ( 2 !)* MI> - k + 1 ) 


(23) 


for the modified Bessel functions I in (13), for large arguments x, valid in the range — 37 T/2 < 


arg(x) < 7 t/2 . In our limit, the expansion (23) turns into the asymptotic approximation 

1 


I,Ax) ~ 


V2 


TTX 


exp i x 


— + e xp 


-x 


1 \ . V 2 \ 

v H— \ m -|- 

2 ) 2 xj 


(24) 
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FIG. 4. (Adapted from [26]) Coordinate system in which the magnetic string, carrying the flux 
<h, pierces the origin of the two-dimensional plane, and x! is aligned along the positive x-axis. For 
processes from x' to x". small angles around ip^ = 0 define the backward direction, and the same 
applies for tpf = 0 for processes in forward direction. 


as we will see in the following. For this purpose, we rewrite the sum in the exact propaga¬ 


tor (13) as 


f h mr"r' 1 _ ^_ 

2 ^ \i+x/m ( 7 t "-t'W 1 e " ^ 


l=—oo 




(25) 


where 


(fib ■— (p" — <p') mod27r G (—7T, 7r), (26) 

see figure [4], such that eventual excess integer multiples of 2n in the angular difference 


appearing in the exponent in (25) are evaluated beforehand. The intention is now to send 


h , as it appears in the index and in the denominator of the argument of the modified Bessel 


functions in (25), to zero in the asymptotic expansion (23) with v = |/ + y/h| and x = —iz 


as defined in (15), except for the exponentials in front of the infinite sums therein. Also in 


the numerator of the exponential of (25) we let h go to zero. The prefactor 1/^/2tix in (23) 
turns into 

1 


,— = \/j— (i + 1 ) - M) , 
V^x V h TTK 


(27) 


where we have introduced the abreviations 


k = 2 mr"r' and t = t" — t'. 


(28) 
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Thus, following the above steps, the product 


P<p = hi (29) 

turns into a continuous variable, and the sum over the quantum mechanically allowed values 


l of the canonical angular momentum in (25) turns into an integral over p^ as given by (29), 

dp^. 


h (/ +1) - hi -> / 

J — c 


(30) 


l=—oo 


In order to see now how (23) turns into (24), we write out the first few terms of the first 


sum in (23) 


£ 


(-1)* T{y + k+\) 1 (-l) c 


^ ( 2x) k k\Y{y - k+\) 0! (2a;) 0 

, 1 (- 1) 1 


1! (2a;) 1 
1 (- 1) 2 
2! (2a;) 2 
1 (- 1) 3 
3! (2a;) 3 


9 

16 

225 259 


64 


+ 


16 



V 

1 

~4 

+v 

5 2 

— V 

2 

+v 

35 4 
4 V 

+1/ 


(31) 


which continues in this manner. Except for the alternating sign, the second series in (23) 


is of the same structure. Since the last column in the brackets in (31) coincides with the 


power series of an exponential, we find that (23) reduces to (24) if all contributions to (31) 


but this last column were to vanish. Upon insertion of the argument x = n/2ith, and of the 


index v — \l + x/h\ of the modified Bessel functions in (25) into (31), we find 

(p v + x)° 


i t 
h hi 
i t 
h hi 
i t 
h n 
i t 
h n 


225 h 6 
64 


+ 


9 h 4 

TeF 

259 h A 


16 


(/V + xf 


_h? 

4 
5 h 2 
2~ 

35 h 2 


+(pv + x¥ 


(Pv + x ) 2 +(Pv + x) 4 ^J 
(P<p + X) 4 +(P V + X) 6 ) + • • ■, 


( 32 ) 
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meaning that the columns prior to the last one are suppressed with respect to the latter 
at least by the order h 2 , and indeed vanish in the limit which we discussed in connection 


with (29) and (30). Thus the contribution to (13) from the first series in (23) turns in this 
limit into 


K _ m 

x ~ 2mh (t" - t') 


exp 


x 


i m (r" — r') 2 


h 2 t" - tj 


T— [°° d Pip exp 


'h ( ^ (?v + x) 


(33) 


The two parts of (23) essentially differ by the extra factor exp (— (v + iti) which comes 

with the second sum. Due to this factor the contribution of the second series in (l23l) vanishes 


in our limit: When splitting the sum in (25) at l x , determined by (17), one obtains two sums 
where in each the index of the modified Bessel functions enters through the extra factor which 


comes with the second series of (23), and contributes a phase of ±7 rl. Except for their sign, 
the resulting contributions of even and odd l coincide in the limit as outlined above, and 
cancel. 


After quadratic completion and evaluation of the Fresnel integral in (|33|) we obtain 

2 


K _ m 

x ~ 2 ?Tin (■t" - t') 


exp 



m (r" — r 
~2 t" -1' 


H- r"r' ^ 


m 

~2 


t" - t' 


XVb 


(34) 


Close to the backward direction, i.e. for small values of (fib, insertion of (26) into (20), fol¬ 


lowed by a truncation of the cosine therein after the second order, produces the phase of (34). 


Therefore (34) and (22) coincide in the region of validity of this expansion, where our limiting 
procedure thus directly leads to the semiclassical approximation of Van Vleck/Gutzwiller. 


Since the integral in (33) is quadratic in p^, also an evaluation by the method of stationary 


phase produces the exact result (34). The stationary point is in this case given by 

<Pb 


p^ = mr"r' 


t" -1' 


X 


(35) 


which, close to the backward direction, coincides with (21) such that our method here 


shows how the main contribution to the integral in (33) results from those canonical angular 
momenta which are close to the classical value. 


In order to determine the asymptotic behavior of (13) close to the forward direction we 
repeat the above steps, but only after insertion of 


iff = ±1t - tpt, € (- 7 T,TT), 


( 36 ) 
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see figure [4j where the sign coincides with the sign of ipb, into (25). Due to the extra term 


±7r in the definition (36), which enters the exponential in (25) as a phase, the same limiting 


procedure as in the above derivation of (34) now suppresses the first series in (23), while the 
second series renders 


m 


K ~ - 

x 2 t rih{t" - t') 


exp 


x 


i m (r" + r'Y 
h~2 I" - V 


kh [L x dPr exp (i U iPr+xf “ 

d'P v ex P ( f x ( “ (Pv + xf - PvPf + XV 


—x 


(37) 


The hnite integration bounds in (37) are determined by multiplication of the constraints on 


L stated in (17) by h, and sending the latter to zero. Since the kinetic angular momentum is 


related to the canonical angular momentum by + X, in the above limit all positive 


kinetic angular momenta contribute to the first integral in (37), whereas the negative kinetic 
angular momenta contribute to the second. 


Due to the hnite integration bounds, evaluation of (37) via the method of stationary 


phase in this case constitutes an approximation. This produces 


K _ m 

x ~ 2mh (t" - if) 


exp 



m (r" + rf m ... p 2 f . . 

-r r ---v (<p f — 7 t) 

2 t" - tJ 2 t" -1' } 


(38) 


for (ff > 0, and 


m 


K v ~ -—-- exp 

x 2 mh (t" - t') y 



77i (r" + r'f rn„,p 2 f 

- r r ---h v {(pf + n) 

2 t" -t' 2 t" -1' ’ 


(39) 


for iff < 0. Close to the forward direction, both (38) and (39) coincide with (22) in the 
respective angular range, such that also in this case the connection to the semiclassical 
propagator of Gutzwiller is established. After the stationary phase approximation leading 


to (38) and (39), respectively, the continuity of (37) is thus no longer resolved. The stationary 


point is in both cases given by 


u / Pf 

Pu, = mr r ---v 

w t „ _ t , v 


(40) 


which also here, close to the forward direction, coincides with the classical value (21) of the 
canonical angular momentum. 
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To discuss in section IV the interference process which smoothes out the discontinuity 
between ( 38p9 ) close to the forward direction, we render the integrals in (37) dimensionless 
and perform a quadratic completion to cast the latter into the final form 


K 


m I i I m ir" + r') m .. . <pi 

~ -exp | — | —-—--- —r r 


x 2n iftt 


h\ 2 t"-t' 


2 t" - 1’ 

1 


x [ exp ( --x (vr - ip f ) ) I dx exp ( ix 2 ) 

^ \ ITT J — ^z/2ifj 

, i . A 1 

+ exp +~x [7T + ipf) —j= / dx exp [ix ) 

ft J V in J -oo 


(41) 


Indeed, from our observation following (37) we find that the contributions of the posi¬ 


tive/negative kinetic angular momenta acquire Dirac’s magnetic phase as if encircling the 
solenoid not more than once in a counterclockwise/clockwise rotational sense. 

Note that nowhere in our derivation of ( 34|41 ) did we explicitly require that tpb or iff be 
close to zero. This requirement must be intricately related to our evaluation of the extra 
phases ±7rZ, since this determines the range of validity of the expressions which result after 
ft is sent to zero. 


IV. INTERFERENCE PROCESS IN THE FORWARD DIRECTION 


In a coordinate system in which x' coincides with the x-axis as in figure |4| the free 


propagator (18) can be split into two half-waves which on their way to x" in the forward 


direction pass the solenoid along the y-axis in an opposite rotational sense ffl- Thereby, 
the half-wave which passes the solenoid counterclockwise via the positive y -axis is given 
by 


Kq =K 0 


J — y/2zsin((ff/2) 


dx exp (ix 2 ) , 


(42) 


see figure 5a The half-wave passing the solenoid clockwise via the negative y -axis is related 


to (42) by K 0 (ip f ) = Jl 0 + (- 99 /), see figure pbl and their sum Kq + K 0 yields the free prop 


agator (18). The additional ad-hoc assumption that these contributions acquire magnetic 
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FIG. 5. Real part of the half-waves K q~ and Kq bending around the origin in a counterclockwise (a) 
and clockwise (b) rotational sense, respectively, divided by the free propagator factor Kq , where x = 
zcosifb, y = zsunpi,. Amended accordingly by Dirac’s magnetic phase, the superposition of these 
half-waves describes the forward Aharonov-Bohm interference process ( 4lf43 ) in the asymptotic 
limit. The sum Kq + Kq amounts to the free-particle propagator I\q. 


phases as in (J22|) then leads to the expression [30] [36] 

K x ~ . 1 


K 0 [ exp ( --y(vr - <p f ) 


+ exp ( + l -x (tt + <p f )j —U J 


J — y/2z sm((ff/2) 

— y/2z sm(cp f /2) 

dx exp ( ix 2 ) 


dx exp (ix 2 ) 


(43) 


which was originally proposed for the forward direction, where it coincides with (41) from 


our novel derivation. The approximate propagator (43) also appears as the lowest order 


of an asymptotic series expansion of (13) |38], where it is shown that the accuracy of (43) 


increases in the entire angular range as z increases. (43), however, exhibits a discontinuity 
as iff —> ±7r, i.e. in backward direction, that vanishes only as z —> oo. The discontinuity 
is, nonetheless, of the same order as the error of the approximation, and whether or not 


there exists an asymptotic approximation of (13) based on (42) that is smooth in the entire 
angular range and containing only the contributions of maximally two topologically distinct 


paths is to date unknown [.39]. For half-integer flux (in units of h), (43) coincides with the 


exact propagator (19). We will elaborate on this point in the next section. 
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(a) z = 7r 





FIG. 6. Comparison of the normalized probability amplitude \K x | 2 /| Kq\ 2 of the exact propaga¬ 


tor (16) (blue, full line) to that of the split-wave approximation (41) (yellow, dot-dashed), centered 
around the forward direction <pf = 0, for y/h = 0.25. The agreement improves as £ increases (a-c). 
Analogous observations hold for the normalized phase arg(/\ v / Kq) (d). 


In figures [6] and [Tj we compare our forward (41) and backward (34) approximations, 
respectively, to the integral representation (16) of the exact propagator (Jl3|. Figure [6] shows 


a comparison of the modulus squared of the full propagator (16) to the semiclassical forward 


approximation (41), both divided by the modulus squared of the free-particle propagator, 
\K X \ J /|A'o| 2 - The expression yields the enhancement of the probability density to find a 
particle prepared at x' at position x" with respect to that of the free evolution, and depends 


only on three parameters: z, as defined in (15), increases with the distance of the preparation 


point x 1 and of the detection point x" from the solenoid, and decreases with increasing 
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—7T -7T/2 0 7T/2 7T 



FIG. 7. Comparison of the normalized probability amplitude |A7J 2 /|A”o| 2 of the exact propa¬ 


gator (16) (blue, full line) to that of the semiclassical (backward) approximation ( 22/34]) (yellow, 
dot-dashed), centered around the backward direction pb = 0, for x/A = 0.25. Comparable accuracy 


to the forward approximation (41) is achieved only for larger z-values (c). 


propagation time, (fif measures the angular distance from the forward direction, x quantifies 
the flux. 


Semiclassical and exact expression show good agreement around the forward direction 


cpf = 0 already for rather small values of z, and the accuracy of (41) increases with z, see 


figure 6a 6c The same applies to the normalized phase arg(A" x /A' 0 ), see figure 6d Around 
the backward direction <pb = 0, comparable agreement between the exact propagator (16) 
and the semiclassical approximation ( 22p4 ) is achieved only for larger z values (of the order 
10 4 7 t), whereat | K x \ 2 /\ A" 0 | 2 converges to unity, see figure [T} 
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FIG. 8. Whenever (41) (yellow, dot-dashed) constitutes a suitable approximation to (16) (blue, full 


line), I-R'xp/I-K’ol 2 scales universally in the angular part R Vf /h = 2 of Hamilton’s principal 

function in units of h. Increasing the square root of z by an arbitrary amount contracts the 
angular range of the interference pattern by the corresponding factor. Ultimately, the width of 


the interference pattern vanishes as yfz —> oo, giving rise to the discontinuity of (22) in forward 
direction (x/h = 0.25). 


In the closing of this section we point out yet another property of the half-wave approx¬ 


imation (41), which facilitates the visualization of the Aharonov-Bohm propagator (13) in 


the asymptotic limit: The normalized probability density depends only on a product of 


the parameters z and iff, which can be inferred neither from (13) nor from (16). Thus, 


whenever (41) is in good agreement with these expressions, we expect that the associated 


interference pattern scales universally in the parameter z^j2 which, close to the forward 
direction, constitutes the angular contribution to Hamilton’s principal function ([3]) of a free 
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particle, 


„ m ... ip 2 * 
Rip f = — t r —-—, 

2 t"-t r 


( 44 ) 


in units of h (see also ( 38p9 |). Thus, if 0 is multiplied by the square of an arbitrary factor, 
we expect the same interference pattern in an appropriately rescaled angular range. This 
is demonstrated in figure [8] where we compare the prediction of (16) to the scale-invariant 


result (41), for several values of z. 


The above implies that, as the width of the interference pattern decreases with the square 
root of z, and vanishes as yfz —> oo (or, equivalenty, y/h —> 0)|30; (31], 40j, it does so in an 


invariant shape. If is positive, the bounds of the integrals in (41) go to —oo as y/h —* 0, 
such that the first integral yields unity while the second vanishes. In turn, only the second 
integral survives for negative ipf. Therefore, in both cases only that contribution to the 
propagator remains which encircles the solenoid in the same rotational sense as the classical 


trajectory associated with semiclassical approximation (22). The latter is thereby recovered 


also in forward direction, exactly as upon evaluation of the integrals in (37) by the method 
of stationary phase, which thus corresponds to neglecting h against all finite values of R Vf . 


Finally, for iff = 0, (43) and (41) can be evaluated directly and yield, independently of the 
value of z, 


(45) 


as first observed in 0Dt 

Our above considerations demonstrate how phase jumps of the semiclassical approxima¬ 
tion (J2| occur entirely in the shadow region defined by the classical trajectories after direct 
passage through the magnetic string. This is also the case, e.g., for the object wave in the 
geometrical shadow of the solenoid in the electron-holographic experiments of Ham. 
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V. RELATION TO THE WHIRLING-WAVE REPRESENTATION 


In this section we discuss the relation of our asymptotic expressions (34) and (41) to the 
whirling-wave representation first introduced in reference HU, 


K = _—_ 

x 2 t rih (t" - f) 


exp 


i m r" 2 + r' 2 


h 2 t" -tr 


x 


E 


dXI\ 




—i Z ) e i ( x -x/ h )(v"-v'+ 27rn ) 


(46) 


which results from the exact Aharonov-Bohm propagator (13) upon invoking the Poisson 
summation formula 


The n-th whirling wave 

/»oo 

T n (Z, (fi) = I 


dXI\ M (- iz ) e iX ^ +2 ™\ 


(47) 


which corresponds to the n-th contribution to the sum in (46), thereby acquires Dirac’s 


magnetic phase exp(— [}p" — <p' + 2im) /h). This property ensures the single-valuedness 


of (46) for arbitrary angular differences <p" — ip 1 HD- For our comparison, we state here 


various aspects of the whirls (47) which, including their asymptotic behavior, have been 


extensively covered in reference 

In figure [9] we plot several of the whirling waves, for different values of 0 . Each of the 


whirls (47) is centered on a particular instance of the backward direction at even integer 


multiples of 7r (dotted vertical lines in figure [9]), and adjacent whirling waves give the largest 
contribution in the forward direction, at odd integer multiples of tt (dashed vertical lines 
in figure [9]). Thus it is common practice to suggest neglecting all but two [28; EH E3H12] 
or three [43| of the whirling waves, which constitutes a sound approximation of the exact 


Aharonov-Bohm propagator (13), in a specific angular range. Such approximations, however, 


feature a small discontinuity that cannot be avoided if any number of the whirling waves (47) 


are excluded from the sum in (46) [43]. Figures [9a] and |9b] show that for increasing z the tails 
of the whirling waves retract, that is the overall contribution of each whirling wave increases 
in a 27r interval centered on its maximum and decreases outside of this interval [32J. In the 
limit z —y oo, the modulus squared of the whirling waves eventually converges to unity inside 
this interval, returning the value 1/2 at the edges in the forward direction, and vanishing 


outside, see figure 9c This corresponds to the semiclassical approximation (22) amended 


for the value (45) in the forward direction iff = 0. 
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0.5 


0 tt 2n 3n 4 tt bn 6 n if 

(c) 0 = 10 1O 7r 


FIG. 9. Amended for Dirac’s magnetic phase, the whirling waves (47) in (46) give rise to an exact 


representation of the Aharonov-Bohm propagator (13) [41 j. (a,b) Each of the whirls is centered on 
a specific instance of the backward direction (gray dotted), and the tails of the whirls retract as z 


increases [43]. Our semiclassical limit (30) shows that, asymptotically, only one whirl contributes in 


the backward direction, see (34), and two whirls in the forward direction (gray dashed), see (41). 
This means that in this limit each whirl cannot extend beyond the backward direction of the 


respective adjacent whirls, (c) In the limit z —> oo the semiclassical propagator (22) is recovered, 


with the value (45) in the forward direction (not shown). 
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In fact, our semiclassical analysis in section Em above shows that in the backward di¬ 


rection the contribution of only one of the whirling waves remains, cf. (34), and in the 


forward direction the same is true for only two of the whirling waves, cf. (41), in our semi¬ 
classical asymptotic limit 


J. Each whirling wave thus does not contribute beyond the 
backward direction of the respective adjacent whirls in this limit. We therefore surmise 


that if an asymptotic approximation of (13) exists, smooth in the entire angular range and 


containing the contribution of only two of the whirls (47), asymptotically, one of these 
whirls must smoothly vanish as the backward direction is approached, while the other whirl 
must smoothly converge to coincide with the semiclassical approximation in the backward 


direction (34). 


For half-integer flux y/A, even and odd contributions to the sum in (46) contribute with 
the same magnetic phase. Thus we find 


m 


K — - 

x 2mh (t" - t') 


x 


exp 


exp 




- h xW 


exp 




i m r" 2 + r' 2 
h~2 t"-t' 


£j 

i even ^ 00 

' T 2vr)] ^ / 

' nodd J ~ 


d\ I\ 


l-M 


d\ L 


-iz e 


(p'-\- 27 rn) 


|A| 


-iz e 


i A ( (p ff —(f 1 + 27 rn) 


)■ 


(48) 


and indeed the first line of (48), together with the sum in the second line, yield the upper 


half-wave Aq defined in (42), while the first and third line, accordingly, yield the lower 


half-wave K 0 m- Thereby, the connection to ( |19| ) is established, which shows why the 
exact half-integer flux propagator already coincides with its semiclassically asymptotic form 


in the forward direction, compare (41), for arbitrary z- values. 


VI. CONCLUSIONS 

In order to address discontinuities which arise when the semiclassical approximation of 
Van Vleck/Gutzwiller is applied to approximately describe magnetic strings, we have in¬ 
troduced a novel semiclassical limit for the Aharonov-Bohm propagator. In this limit, the 
discrete sum over the quantum mechanical canonical angular momenta turns into integrals 
over a continuous variable, and the stationary points of those integrals coincide with the 
classical value of the canonical angular momentum. While our limit directly leads to the 














23 


semiclassical approximation of Van Vleck/Gutzwiller in the backward direction, it generates 
in the forward direction an asymptotic split-wave expression of two contributions passing the 
string on either side. Thereby, the positive/negative kinetic angular momenta acquire Dirac’s 
magnetic phase as if passing the string not more than once in a counterclockwise/clockwise 
rotational sense. We pointed out that the interference pattern generated by our split-wave 
approximation scales universally in the angular part of Hamilton’s principal function in units 
of Planck’s constant. Divergence of this ratio gives rise to the above-mentioned disconti¬ 
nuities, which corresponds to an evaluation of the canonical angular momentum integrals 
by means of the stationary phase approximation, and eventually leads to the expression of 
Van Vleck/Gutzwiller also in the forward direction. Our approach hints at features of a - 
hypothetical - globally smooth asymptotic approximation, whence we surmise that the form 
of such a function coincides in the respective regions with our asymptotic expressions for 
the forward and backward directions. We have shown that of a representation of the exact 
propagator in terms of whirling waves remain only one contribution in the backward direc¬ 
tion and two contributions in the forward direction, in our limit. In the case of half integer 
flux in which the propagator can be written in a closed form, every other whirl contributes 
with the same magnetic phase factor, directly giving rise to one of the half-waves which are 
for an arbitrary value of the flux encountered in the forward direction only in the asymptotic 
limit. 
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